LIFE EXTINCTIONS DUE TO 
NEUTRON STAR MERGERS 

Arnon Dar, Ari Laor, and Nir J. Shaviv 

Department of Physics and Space Research Institute 
Israel Institute of Technology, Haifa 32000, Israel. 

'<6 '. 

, Abstract 

o\ ■ 

y—i ' Cosmic ray bursts (CRBs) from mergers or accretion induced collapse 
^ ■ of neutron stars that hit an Earth-like planet closer than ~ 1 kpc from 
O , the explosion produce lethal fluxes of atmospheric muons at ground level, 
I underground and underwater. These CRBs also destroy the ozone layer 
' and radioactivate the environment. The mean rate of such life devastat- 
es! : ing CRBs is one in 100 million years {Myr), consistent with the observed 
\ 5 "great" extinctions in the past 600 Myr. Unlike the previously sug- 
^> ■ gested extraterrestrial extinction mechanisms the CRBs explain massive 
. life extinction on the ground, underground and underwater and the higher 
\ survival levels of radiation resistant species and of terrain sheltered species. 
' More distant mergers can cause smaller extinctions. Biological mutations 
^ ■ caused by ionizing radiation produced by the CRB may explain a fast 
Q>^ i appearance of new species after mass extinctions. The CRB extinction 
' predicts detectable enrichment of rock layers which formed during the ex- 
p i' tinction periods with cosmogenically produced radioactive nucleides such 

o 



as ^^^I, ^^^Sm, 205p|j and ^'^^Pu. Tracks of high energy particles in 
|H \ rock layers on Earth and on the moon may also contain records of intense 

^ ■ CRBs. An early warning of future extinctions due to neutron star mergers 

can be obtained by identifying, mapping and timing all the nearby binary 
neutron stars systems. A final warning of an approaching CRB from a 
^ ' nearby neutron stars merger will be provided by a gamma ray burst a few 

. days before the arrival of the CRB. 
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1 Introduction 



The early history of single-celled organisms during the Precambrian, 4560 to 570 
million years {Myr) ago, is poorly known. Since the end of the Precambrian the 
diversity of both marine and continental life increased exponentially. Analysis of 
fossil record of microbes, algae, fungi, protists, plants and animals shows that this 
diversification was interrupted by five major mass extinction "events" and some 
smaller extinction peaks (1) . The "big five" mass extinctions occurred in the 
Late Ordovician, Late Devonian Late Permian, Late Triassic and end-Cretaceous 
and included both marine and continental life. The largest extinction occurred 
about 251 Myr ago at the end of the Permian period. The global species ex- 
tinction ranged then between 80% to 95%, much more than, for instance, the 
end-Ordovician extinction 439 Myr ago which eliminated 57% of marine gen- 
era, or the Cretaceous- Tertiary extinction 64 Myr ago which killed the dinosaurs 
and claimed 47% of existing genera (2). In spite of intensive studies it is still 
not known what caused the mass extinctions, how quick were they and whether 
they were subject to regional variations. Many extinction mechanisms have been 
proposed but no single mechanism seems to provide a satisfactory explanation 
of both the marine and continental extinction levels, the biological extinction 
patterns and the repetition rate of the mass extinctions (1,2). These include 
astrophysical extinction mechanisms, such as meteoritic impact that predicts the 
iridiun anomaly which was found at the Cretaceous/Tertiary boundary (3) but 
has not been found in all the other extinctions (4), supernova explosions (5) and 
gamma ray bursts (6) which do not occur close enough at sufficiently high rate to 
explain the observed rate of mass extinctions. In this paper we propose that cos- 
mic ray bursts (CRBs) from neutron star (NS) mergers at distances closer than 
1 kpc from planet Earth produced lethal fiuxes of atmospheric muons at ground 
level, underground and underwater, destroyed the ozone layer radioactivated the 
atmosphere and the surface of the planet and polluted it with radioactive iso- 
topes. Such CRBs can explain the large mass extinctions on planet Earth in 
the past 600 Myr. Unlike other extraterrestrial extinction mechanisms that have 
been proposed before, CRBs are much more effective in extincting life and can 
explain massive extinction underwater and underground, the higher survival lev- 
els of radiation resistant species and terrain sheltered species, and the observed 
rate of the large mass extinctions on Earth. More distant galactic mergers may 
have caused some of the smaller extinctions on our planet. Biological mutations 
induced by ionizing radiation could have caused a fast appearance of new species 
after mass extinctions. 
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2 CRBs From NS Mergers 



Three NS-NS binaries are presently known within the galactic disc; PSRs B1913+16 
(7) at a distance of D ~ 7.3 kpc, B203+46 (8) at D ~ 2.3 /cpc, and B1534+12 
(9) at D ~ 0.5 kpc, and PSR 2127+llC (10) in the globular cluster M15 at 
D ~ 10.6 kpc. Although unseen as radio pulsars, the companion stars in these 
systems have been identified as neutron stars from their mass which is inferred via 
measurements of the relativistic periastron advance using standard pulse timing 
techniques (11). These binary neutron stars loose energy by gravitational wave 
emission. The continuous energy loss brings them closer and closer until they 
finally merge. Their merger times due to gravitational radiation have been cal- 
culated (12) from the observed binary period, orbital eccentricity and the masses 
of the pulsar and its companion. They arc 3.01 x 10*^ yr, t ^ tn ■ 2.73 x 10^ yr 
and 2.20 x 10® yr, respectively, where tn = l/H is the Hubble age of the Uni- 
verse. It has been estimated that there are ~ 240 potentially observable NS-NS 
binaries in the MW disc and their merging rate is (13) Rmw ~ 2 — 3 x {Myr)~^. 
Such NS-NS mergers release enormous amount of gravitational binding energy, 
~ MqC^, in few ms in the form of gravitational waves, neutrinos and kinetic 
energy of ejected material (14). The ejecta probably consists of high Z material 
because the crust of neutron stars consists of iron-like nuclei and because even 
if the ejecta is initially made of neutrons or dissociated nuclear matter, upon 
expansion nucleosynthesis transforms it quickly into the most bound, iron-like 
nuclei (15). The kinetic energy of the ejecta can be estimated from numerical 
simulations of NS-NS mergers (16). Such estimates are highly questionable since 
numerical simulations have difficulties in reproducing even Type II supernovae 
explosions, let alone the much more complicated NS-NS mergers. However, a 
rough estimate of this energy can be obtained directly from observations. Typi- 
cally in SNell, IOMq are accelerated to a final velocity of (12) v ~ 10^ km s~^, 
i.e., to a total final momentum P ~ IQMqV. Since core collapse is not affected 
directly by the stellar envelope, we assume that a similar impulse, / Fdt ~ P, 
is imparted to the ejected mass in merger/ AIC collapse of NS (this may not 
be the case if energy and momentum deposition of neutrinos in the envelope is 
driving the explosion). If the ejected mass is much smaller than a solar mass, 
AM <C Mq, it is accelerated to a highly relativistic velocity. Its kinetic energy 
then is, Ek ~ Pc ~ 5 x 10^^ dil erg' in a solid angle dQ, which is injected into the 
interstellar medium as a very powerful cosmic ray burst (CRB). 

3 Indirect Evidence From GRBs 

It was shown (17) that highly relativistic ejecta from NS-NS mergers can boost 
and beam star light in dense stellar regions into cosmological gamma ray bursts 
(GRBs) very similar to those observed (18), provided its kinetic energy is ~ 
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5 X lO^^dn erg and T = E/mc^ ~ 1000, where dQ > l/AF^ is the sohd angle of 
the beamed emission. In particular, such CRBs in dense stellar regions explain 
remarkably well (17) the burst size, bimodal duration distribution, complex light 
curves and spectral evolution of GRBs. Moreover, the estimated rate of NS-NS 
mergers in the Universe (19) is consistent with the observed rate (18) of GRBs and 
the cosmic ray injection rate in the MW that is required to maintain a constant 
cosmic ray flux in the MW: 

The escape rate of cosmic rays from the MW requires an average injection 
rate of Qcr ~ 5 x 10^° erg s~^, in high energy cosmic rays above GcV, in order 
to maintain a constant energy density of cosmic rays in the MW (20), where we 
used a total mass of gas ~ 4.8 x IO^Mq in the MW disc deduced from cosmic 
ray production of galactic 7 rays. This injection rate is consistent with that due 
to NS-NS mergers in the disc of the MW, 

Qc^ ~ RmwEk ~ 5 X 10^° erg s''. (1) 

4 Attenuation of CRBs 

The typical burst size and duration of GRBs suggest that the energy release by 
NS-NS merger in high energy cosmic rays into a solid angle dQ > l/4r^ in the 
direction of Earth is ~ 5 x 10^'^d^ erg and that T = E/mc^ ~ 1000 (17). Such 
CRBs cannot be attenuated significantly by the column densities of radiation and 
gas in the galactic disk. This is unhke the non relativistic ejecta in supernova 
(SN) explosions, which is attenuated by the interstellar gas over a distance of a 
few pc by Goulomb collisions: Moderately energetic charged particles, other than 
electrons, lose energy in neutral interstellar gas primarily by ionization. The 
mean rate of energy loss (or stopping power) is given by the Bethe-Bloch formula 

- dE/dx « (47rZ2ne/?-2eVmec2)[(l/2)/n(2mec2/?2f r„„J/2) _ p^]^ (2) 

where Ze is the charge of the energetic particle of mass Mj, velocity /3c and total 
energy E = jMiC^, is the number of electrons per unit volume in the medium 
in atoms with ionization potential /, and T^ax = 2meC^/3^7^/(l + 2'jme/Mi). 
Hence, for SN ejecta with Z ~ 1 and j3 = v/c ^ 10^^ in a neutral interstel- 
lar medium with a typical density of uh ~ 1 cm^^ the stopping distance of 
the ejecta due to Coulomb interactions is x ~ E/2{dE/dx) ~ 0.7 pc. If the 
interstellar gas around the SN is ionized, then I has to be replaced by e^/Ro 
where Rd = [kT / ATie^ne{Z + 1))^^^ is the Debye screening length. Thus for an 
ionized interstellar gas with kT ~ 1 eV the range is a; ~ E/2{dE/dx) ~ 6 
if t; ~ 10000 km s~^. In fact, the stopping of the ejecta by Goulomb interac- 
tions is consistent with observations of SN remnants, like SN 1006 (21), while 
the assumption that the ionized interstellar medium is glued to the swept up 
magnetic field seems to be contradicted by some recent observations (22). The 
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stopping power of the interstellar medium for highly relativistic iron-like nuclei, 
which dominate the CRB composition, is due mainly to nuclear collisions and 
radiative energy losses and not to ionization losses. The range of iron-like nuclei 
with r ~ 1000 in hydrogen is approximately ~ lO^^cm in a typical interstellar 
density of hh ~ 1 cm~^, while their range due only to ionization energy loss is 
~ 10^^ cm. 

Although the galactic magnetic field, ~ 3 — 5 x 10~^ Gauss, results in a 
Larmor radius of = j3Vmc/qH ~ 10^'^ cm for protons with V = 1000, it cannot 
deflect and disperse the CRB particles from a NS-NS merger at distance smaller 
than ~ 1 kpc from the explosion. This is similar to the case of highly relativistic 
jets from active galactic nuclei which reach distances of hundreds of kpc without 
significant defiection. Because of its high energy density, the CRB sweeps up 
the interstellar magnetic field as long as its distance from the explosion is smaller 
than D ~ {<oEk/H^Y/^ ~ 2 kpc, where its energy density becomes comparable to 
that of the swept up magnetic field. (The energy densities of the galactic cosmic 
rays and the magnetic field in the galactic disc are e^^ ~ ~ 10~^^ erg cm~^, 
respectively). Thus, CRBs expand ballistically, like jets from active galactic 
nuclei, up to a distance D ~ 1 kpc. 

5 Mass Extinction By CRBs 

We assume that the ambient interstellar gas is not swept up with the CRB. If 
it were, then the CRB would not expand beyond ~ 10 pc. In SN explosions 
collective modes are invoked as the source of the couphng of the ejecta to the 
interstellar medium, required in order to attenuate the SN debris. As mentioned 
above, binary Coulomb interactions arc sufficient to produce the observed cou- 
pling in SN explosions, and coupling through collective modes is not necessarily 
present. 

Unattenuated CRBs from NS-NS mergers can be devastating to life on nearby 
planets. At a distance of 1 kpc the duration of unattenuated CRB is 

St ~ D/2cr^ r^lday-2 months (3) 

for typical values of T between 1000 and 100, respectively. The time integrated 
energy fiux of the CRB at D ~ 1 kpc is, typically, 3 x 10^^ GeV cm'"^. Thus, the 
energy deposition in the atmosphere by the CRB is equivalent to the total energy 
deposition of galactic cosmic rays in the atmosphere over ~ 10^ yr. However, the 
typical energy of the cosmic rays in the CRB is ~ 1 TeV per nucleon, compared 
with ~ 1 GeV per nucleon for ordinary cosmic ray nuclei. Collisions of such CRB 
particles in the atmosphere of a planet generate atmospheric cascades where a 
significant fraction of the CRB energy is converted into "atmospheric muons" 
through leptonic decay modes of the produced mesons. Most of these muons do 
not decay in the atmosphere because of their high energy, unlike most of the 
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atmospheric muons that are produced by ordinary cosmic rays . The average 
number of high energy muons produced by nucleons of primary energy Ep, which 
do not decay in the atmosphere and reach sea level with energy > at zenith 
angle 6 < 7r/2, is given approximately by (23) 

< {O.OU5E^[TeV]){Eo/E^f\l - E^/Eof/cose . (4) 

Thus a CRB with energy of about 1 TeV per nucleon at a distance of 1 kpc 
produces at sea level a flux of atmospheric muons of 

7^(> 3 GeV) ~ 10^^ cm-^. (5) 

Such muons deposit energy in matter via ionization. Their energy deposition 
rate is (24) —dE/dx > 2 MeV g~^cm~^. The whole-body lethal dose from 
penetrating ionizing radiation resulting in 50% mortality of human beings in 
30 days (24) is < 300 rad ^ 2 x W°/{dE/dx) - 10^° cm-^ where dE/dx 
is in rate is in MeV g~^cm~^ units. The lethal dosages for other vertebrates 
can be a few times larger while for insects they can be as much as a factor 20 
larger. Hence, a CRB at D ~ 1 kpc which is not significantly dispersed by 
the galactic magnetic field produces a highly lethal burst of atmospheric muons. 
Because of muon penetration, the large muon flux is lethal for most species even 
deep (hundreds of meters) underwater and underground, if the cosmic rays arrive 
from well above the horizon. Thus, unhke the other suggested extraterrestrial 
extinction mechanisms, a CRB which produces a lethal burst of atmospheric 
muons can explain also massive extinction deep underwater and why extinction 
was higher in shallow waters. 

Although half of the planet is in the shade of the CRB, planet rotation exposes 
a larger fraction of the planet surface to the CRB. Additional effects increase the 
lethality of the CRB over the whole planet. They include: (a) The pollution of 
the environment by radioactive nuclei, produced by spallation of atmospheric and 
surface nuclei by shower particles. Using the analytical methods of ref. 25, we 
estimate that for an Earth-like atmosphere, the flux of energetic nucleons which 
reaches the surface is also considerable, 

Ip{> 100 MeV) ~ In{> 100 MeV) ~ 10^° cm'^ (6) 

Global winds spread radioactive gases in a relatively short time over the whole 
planet. 

(b) Depletion of stratospheric ozone by the reaction of ozone with nitric oxide, 
generated by the cosmic ray produced electrons in the atmosphere (massive de- 
struction of stratospheric ozone has been observed during large solar flares which 
produced energetic protons (26)). 

(c) Extensive damage to the food chain by radioactive pollution and massive 
extinction of vegetation and living organisms by ionizing radiations (the lethal 
radiation dosages for trees and plants arc slightly higher than those for animals 
but still less than the flux given by eq. 5 for all except the most resilient species) . 
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6 Signatures of CRB Extinction 



The biological extinction pattern: The biological extinction pattern due to a 
CRB depends on the exposure and the vulnerability of the different species to the 
primary and secondary effects of the CRB. The exposure of the living organisms 
to the muon burst depends on the size and duration of the CRB, on its direction 
relative to the rotation axis of Earth (Earth shadowing), on the local sheltering 
provided by terrain (canyons, mountain shades) and by underwater and under- 
ground habitats, and on the risk sensing/assessment and mobility of the various 
species. The lethahty of the CRB depends as well on the vulnerability of the 
various living species and vegetation to the primary ionizing radiation, to the 
drastic changes in the environment (e.g., radioactive pollution and destruction of 
the ozone layer) and to the massive damage and radioactive poisoning of the food 
chain. Although the exact biological signature may be quite complicated, and 
somewhat obscured in fossil records (due to poor or limited sampling, deteriora- 
tion of the rocks with time and dating and interpretation uncertainties because of 
bioturbational smearing) it may show the general pattern expected from a CRB 
extinction. Indeed, a first examination of the fossil records suggest that there 
is a clear correlation between the extinction pattern of different species, their 
vulnerability to ionizing radiation and the sheltering provided by their habitats 
and the environment they live in. For instance, insects which are less vulnerable 
to radiation, were extinct only in the greatest extinction - the end-Permian ex- 
tinction 251 Myr ago. Even then only 8 out of 27 orders were extinct compared 
with a global species extinction that ranged between 80% to 95% (4). Also plants 
which are less vulnerable to ionizing radiation suffered lower level of extinction. 
Terrain, underground and underwater sheltering against a complete extinction 
on land and in deep waters may explain why certain families on land and in deep 
waters were not extinct even in the great extinctions, while most of the species 
in shallow waters and on the surface were extinct (4). Mountain shadowing, 
canyons, caves, underground habitats, deep underwater habitats and high mobil- 
ity may also explain why many species like crocodiles, turtles, frogs, (and most 
freshwater vertebrates), snakes, deep sea organisms and birds were little affected 
in the Cretaceous/Tertiary (K/T) boundary extinction which claimed the life 
of the big dinosaurs and pterosaurs. In particular, fresh underground waters in 
rivers and lakes are less polluted with radioisotopes and poisons produced by the 
CRB than sea waters and may explain the survival of freshwater amphibians. 

Geological Signatures: The terrestrial deposition of the primary CRB nu- 
clides or the production of stable nuclides in the atmosphere or in the surface 
by the CRB is too small to be detectable. In particular, the proposed mech- 
anism cannot explain the surface enrichment at the K/T boundary by about 
3 X 10^ tons of iridium. Alvarez et al. (3) suggested that the impact of ex- 
traterrestrial asteroid, with an Ir abundance similar to that observed in early 
solar system Chondritic meteorites whose Ir abundance is larger than crustal Ir 
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abundance by ~ 10 , could cause the Ir anomaly and explain the mass extinc- 
tion at the K/T boundary. However, no significant Ir enrichment was foTind in 
all other mass extinctions. Moreover, other isotopic anomalies due to meteoritic 
origin have not been found around the K/T boundary; in particular the As/Ir 
and Sb/Ir ratios are three orders of magnitude greater than chondritic values 
but are in accord with a mantle origin (27). Extensive iridium measurements 
showed that the anomaly does not appear as a single spike in the record, indica- 
tive of an instantaneous event, but rather occur over a measurable time interval 
of 10 to 100 kyr or possibly longer (28). That and the high abundance of Ir 
in eruptive magma led to the suggestion that the iridium anomaly is due to a 
global volcanic activity over 10 to 100 kyr at the K/T boundary (28,29) which 
also caused the K/T extinction: Eruptions have a variety of short term effects, 
including cooling from both dust and sulfates ejected into the stratosphere, acid 
rain, wildfires, release of poisonous elements and increase in ultraviolet radiation 
from ozone-layer depletion. But, examination of major volcanic eruptions in the 
past 100 Myr have shown that none of them greatly affected the diversity of 
regional and global life on land or in the oceans (4). 

A CRB will enhance the abundance of stable cosmogenic isotopes in the ge- 
ological layer corresponding to the CRB event, but, the enrichment may be neg- 
ligible compared to their accumulation through long terrestrial exposure of the 
geological layers to galactic cosmic rays prior to the CRB. However, CRB enrich- 
ment of sediments with unstable radioisotopes of mean lifetimes much shorter 
than the age of the solar system, r <^ ^ 4570Myr, but comparable to the 
extinction times, may be detectable through low traces mass spectrometry. In 
particular, fission of long lived terrestrial nuclei, such as ^^^U and ^^^Th, by shower 
particles, and capture of shower particles by such nuclei, may lead to terrestrial 
production of, e.g., ^^^I with r — 15 Myr, ^^^Sm with r = 146 Myr, 205p|-j ^j^j^ 
r = 43 Myr and ^^^Pu with r = 118 Myr, respectively. These radioisotopes may 
have been buried in underwater sediments and underground rocks which were 
protected from further exposure to cosmic rays. The main background to such 
a CRB signature is the continuous deposition by cosmic rays and by meteoritic 
impacts on land and sea. Cosmic rays may include these trace radioisotopes due 
to nearby sources (e.g., supernova explosions) and because of spallation of stable 
cosmic ray nuclei in collisions with interstellar gas. Meteorites may include these 
trace elements due to a long exposure in space to cosmic rays. 

Finally, large enhancement of TeV cosmic ray tracks in magma from volcanic 
eruptions coincident with extinctions may also provide fingerprints for CRB ex- 
tinctions. 
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7 Rate of Mass Extinctions 



Assuming that the spatial distribution of NS binaries in the MW follow the 
distribution of single pulsars (11), 

dN oc e-^'/2i?^g-|.|A^^^^^^ (7) 

with a disc scale length, i?o ~ 4.8 fcpc, and a scale height, h > 0.5 kpc perpendic- 
ular to the disc and independent of disc position, we find that the average rate 
of NS-NS mergers within ~ 1 kpc from planet Earth is ~ 10~^ yr~^. It is con- 
sistent with the 5 big extinctions which have occurred during the last 600 Myr 
in the Paleozoic and Mesozoic eras. The relative strengths of these extinctions 
may reflect mainly different distances from the CRBs. Beyond ~ 1 kpc from the 
explosion the galactic magnetic field begins to disperse the CRB and suppresses 
its lethality. Such CRBs, if not too far, can still cause partial extinctions at a 
higher rate and induce biological mutations which may lead to the appearance of 
new species. 

The galactic rate of SN explosions is ~ 100 y~^. The range of debris from 
SN explosions in the interstellar medium is shorter than 10 pc. The rate of 
SN explosions within a distance of 10 pc from Earth which follows from eq. 7 is 
Rsn{< 10 pc) ~ 10~^%~^. High energy cosmic rays which, perhaps, are produced 
in the SN remnant by shock acceleration, carry only a small fraction of the total 
explosion energy and arrive spread in time due to their diffusive propogation in 
the interstellar magnetic field. Also neutrino and fight emissions similar to those 
observed in SN1987A, at a distance of a few pc cannot cause a mass extinction. 

8 Conclusions 

Cosmic Ray Bursts from neutron star mergers may have caused the massive con- 
tinental and marine life extinctions which interrupted the diversification of life 
on our planet. Their rate is consistent with the observed rate of mass extinctions 
in the past 570 Myr. They may be able to explain the complicated biological 
and geographical extinction patterns. Biological mutations induced by the ion- 
izing radiations which are produced by the CRBs may explain the appearance 
of completely new species after extinctions. A first examination suggests a sig- 
nificant correlation between the biological extinction pattern of different species 
and their exposure and vulnerability to the ionizing radiation produced by a 
CRB. The iridium enrichment around the Cretaceous/Tertiary extinction that 
claimed the life of the dinosaurs and pterosaurs cannot be due to a CRB. It may 
have been caused by intense volcanic eruptions around that extinction. Isotopic 
anomaly signatures of CRB extinctions may be present in the geological layers 
which recorded the extinctions. Elaborate investigations of the effects of CRBs 
from relatively nearby neutron star mergers and their biological, radiological and 
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geological fingerprints are needed before reaching a firm conclusion whether the 
massive extinctions during the long history of planet Earth were caused by CRBs 
from neutron stars mergers. If nearby neutron star mergers are responsible for 
mass extinctions, then an early warning of future extinctions due to neutron star 
mergers can be obtained by identifying, mapping and timing all the nearby binary 
neutron stars systems. A final warning for an approaching CRB from a nearby 
neutron-stars merger will be provided few days before its arrival by a gamma ray 
burst produced by the approaching CRB. 
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